The role of glass dynamics in the anomaly of the dielectric function of solid helium 
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We propose that acousto-optical coupling of the electric field to strain fields around defects in 
disordered *He is causing an increase of the dielectric function with decreasing temperature due 
to the arrested dynamics of defect excitations. A distribution of such low-energy excitations can 
be described within the framework of a glass susceptibility of a small volume fraction inside solid 
■^He. Upon lowering the temperature the relaxation time r{T) of defects diverges and an anomaly 
occurs in the dielectric function e(cj,T) when ujt{T) ~ 1. Since e{u},T) satisfies the Kramers-Kronig 
relation, we predict an accompanying peak in the imaginary part of e{ui, T) at the same temperature, 
where the largest change in the amplitude has been seen at fixed frequency. We also discuss recent 
measurements of the amplitude of the dynamic dielectric function that indicate a low-temperature 
anomaly similar to the one seen in the resonance frequency of the torsional oscillator and shear 
modulus experiments. 

PACS numbers: 67.25.dt, 67.80.B-, 67.80.bd 
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The reported anomalies in resonance frequency and 
dissipation of torsional oscillators [l[ and in the dynamic 
shear modulus [2l at low temperatures are the subject 
of intense studies as they were suggested to be the sig- 
natures of supersolidity. It is now generally agreed that 
disorder plays an essential role for observing these ef- 
fects. The challenge of unambiguously identifying Bose- 
Einstein condensation (EEC) in solid *He arises from the 
fact that crystal defects are required for supersolidity ■ 
These defects exhibit their own dynamics and contribute 
to the observed period of oscillations, shear modulus and 
specific heat in the same temperature and pressure range, 
where supersolidity is expected. Thus any unambiguous 
identification of a supersolid state requires a detailed un- 
derstanding of the behavior of defects. 

To study the role of glassy excitations, which make 
up a small volume fraction of the crystal, we proposed a 
phenomenological framework that captures the dynam- 
ics of defects 0, . It accounts for observed anomalies 
in torsional oscillator, specific heat, and shear modulus 
experiments 0-11] and can be thought of as a distribu- 
tion of crystal defects forming two-level-systems (TLS) 
Q- Possible candidates for the TLS are groups of atoms 
near a crystal defect or pinned segments of dislocation 
lines. These defects can be removed through annealing or 
created by applying shear stress, thus drastically chang- 
ing the mechanical and elastic properties of solid ''He 
Q. A freezing out of these excitations can account for 
the anomalies by postulating a relaxation time diverging 
with decreasing temperature. This is typical of glassy 
systems overcoming an activation barrier which 
results in visco-elastic behavior, long relaxation times, 
hysteresis, as well as effective softening of elastic moduli. 

Given the growing experimental evidence for the role 
of disorder in solid ^He, we ask the question about the 
effects of a glassy component on the dielectric function 
at lowest temperatures. Here we show that glassy dy- 



namics can affect the strain and position of atoms, which 
in turn leads to an additional contribution to the po- 
larization. This contribution causes a decrease of the 
dielectric function at high temperatures. Therefore, we 
propose that the cooperative motion of displaced (out-of- 
equilibrium) atoms near defects, in the glasslike regions 
of solid, creates local strain and thus reduces the polar- 
ization through acousto-optical coupling. Note that this 
effect is not captured by the standard Clausius-Mossotti 
equation for dielectrics, which attributes the change in 
dielectric function e to either a change in polarizability 
a or in mass density p: (e — l)/(e + 2) = (47r/3)(ap/Af), 
with molar mass M. In fact, the glass model yields an ef- 
fect opposite in sign and orders of magnitude bigger than 
a correction in the polarizability, predicted for dipole- 
induced dipole interactions [9| , which was reported a long 
time ago in ''He by Chan et al. (loj . 

The glass model describes the anomaly in e{uj,T) due 
to acousto-optical coupling between strain induced by 
glassy defects and dielectric polarization. It features: (i) 
a decrease of e{uj, T) at high temperatures; (ii) a dissipa- 
tion peak in the same temperature range where the am- 
plitude of e(w, T) changes most significantly. The origin 
for this dynamic behavior is set by the matching condi- 
tion ujt{T) ~ 1 between the frequency uj and relaxation 
time t(T). The key assumption is that a strain field 
in the vicinity of crystalline defects induces changes in 
the polarization, since the magnitude of the local electric 
field sensed by a helium atom is different from the elec- 
tric field at equilibrium atomic position in a crystal free 
of internal stress. 

Very recent measurements of e{uj,T) by Yin and 



coworkers show the increase of dielectric function at 
low temperatures. We illustrate that these data cannot 
be described by the standard Clausius-Mossotti equation 
through a change in mass density or polarizability, e.g., 
due to dipole-induced dipole interactions. Neither the 
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measured change of the mean mass density 5pj p 10~^, 
nor the predicted correction in polarizabihty, which ac- 
tually leads to a decrease of e(w, T) at low temperatures 
0, [13 > can account for the reported change of the di- 
electric function of order 5e/e~ 10~^, when invoking the 
Clausius-Mossotti equation. 

The minimal coupling model between polarization and 
strain field is obtained by expanding the microscopic po- 
larization, p(™icro)^ around its equilibrium value: 

p(micro) (j.) _ p(micro) (-j^^, + (u • V) p("""°) , (1) 

where r = R + u is the atom's position, R is its equi- 
librium position, and u is its displacement. The macro- 
scopic polarization P measured in experiment is related 
to p (micro) ^Yie average over a macroscopic volume 
element v, P = IJ^dr p('"*"°)(r). In the presence of 
a local strain field one finds in linear order of the dis- 
placement P « / dR [P('"'"-°)(R) + (u • V) p(«""'°)(r)]. 
Integration of the first term gives the macroscopic po- 
larization for zero internal strain (a solid in equilib- 
rium), Pq = 47rxoE where xo is the zero-strain sus- 
ceptibility. The second term can be treated through 
integration by parts, that is, / dR (u • V) p(micto) ^ 
Jg dA (u • n) p(™"°) - / dR (V • u) p(micro)^ The u • n 
term vanishes at the surface of v, since the strain field 
is localized within a macroscopic distance. The second 
term, assuming a smoothly varying strain field, modifies 
the polarization by 

SP ^ eu JdR P("""°)(R) = -e,,Po, (2) 

with the dilatory strain = V ■ u (we use the notation 
en = eii -|- 622 + 633). This change in local polarization 
corresponds to a change in dielectric function of 

Sea (w, T) = -Attxo eu (w, T). (3) 

Note that only the diagonal components of the strain 
tensor play a role in this leading order expansion. In 
principle, the shear strain can couple to the electric field 
by considering dipole-induced dipole interactions (Van 
der Waals), which is a higher order efltect. Note that 
our derivation of Eq. ([3]) is equivalent up to leading or- 
der to the Pockels coefficient for acousto-optical cou- 
pling in isotropic or polycrystalline dielectrics, Je^ — 
— [2P44 Cij -I- P12 efefeJy ], where Pki are the reduced 
Pockels coeflticients (of order unity) |12| . 

Next we discuss the origin of the local strain. De- 
pending on the rate of the solidification process of he- 
lium, various amounts and kinds of defects can be frozen 
into the crystal leaving behind relics of the history of 
the phase transition. We postulate that defects in the 
glassy regions of solid form the TLS. The associated 
strain fields are localized near the TLS and exhibit glassy 
dynamics, since nearby atoms reshuffle with an intrin- 
sic time delay. We describe the equation of motion of 



these atoms within linear response theory by including a 
back action term @] . In an isotropic or polycrystalline 
medium the elastic stress tensor = Xijkiduk/dxi 
takes a simplified form with the elastic modulus tensor 
Kjki = >^o6ijSki + po{SikSji + SiiSjk). If the electric field 
couples to local density fluctuations only through dila- 
tory strain, then the important matrix element is the 
Lame parameter Aq of the purely elastic solid. We write 
the displacement to an out-of-equilibrium internal force 
in the presence of the back action as 

p 52 U,{t) + Ao dfu,{t) = /™T(i) + /BA(^)^ 

where Ui is the displacement of an atom in the ith di- 
rection and p is the mass density, f}^"^ is the out-of- 
equilibrium internal force density in the ith direction at 
the defect, is the back action force density that 

describes the time-delayed response of nearby atoms: 

fFHt)= f dt'g{t~t';T) dfu.it'), (5) 

J —00 

where Q is the strength of the back action onto bulk *He. 
Integration of Eq. ^ over the ith direction and Fourier 
transformation to frequency domain yields the dynamic 
strain due to a local dilatory stress an: 

e..(a.,r)-a,,(Ao-e(c.,T))-\ (6) 

We assumed that the back action can be described 
by a distribution P{t) of Debye relaxors, Q{uj,T) — 
dt P{t) [1 — 1/(1 — joj t{T) t)] , with relaxation time 
t(T) of the glassy regions of solid. The specific form 
of t{T) can change quantitatively, but not qualitatively, 
the T-dependence of en. For simplicity, we choose for 
P{t) the Cole-Cole distribution resulting in Q{ijj,T) = 
goAo/ [1 — (iw r)"] . Therefore the corresponding dilatory 
strain of a distribution of glassy TLSs becomes 

en{io,T) = eo(l~.go/[l-(ic^r)"])"\ (7) 

where cq = Cm/Aq at T = and 170 is the glass parameter, 
which depends on the TLS density. From Eq. ([7]) and ([3|) , 
we obtain the change in dielectric function due to local 
strain fiuctuations, 

5en = -4^Xoeo(l-.go/[l-(«wr)"])"'. (8) 

At low temperatures, r — > 00 and en — > cq, hence the 
strain is minimal and the reduction of dielectric function 
due to strain fiuctuations is small. At high temperatures, 
T — > and en — > eo(l — .90)^^ reaches its maximum re- 
sulting in the largest reduction of strain, where solid ^He 
is softest. This is the main result of this study, namely 
that the dielectric function reflects the arrested dynam- 
ics of the glassy components at low temperatures through 
the acousto-optical coupling. Within this approximation 
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FIG. 1: (Color online) Experimental data and calculations 
of the dielectric function vs. temperature. The red circles 
are the ex per imental data of the dielectric function (data by 
Yin et al. The black lines are the calculated amplitude 

and phase lag (dissipation) of e{uj,T). We used parameters 
a = 1.49, eo = 8.88 x 10"*, 5 = 0.21, to = 10.4 ns, A = 1.92 
K, and Tg = -119 mK. 
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FIG. 2: (Color online) Low-temperature pressure deviation 
from lattice contribution of Debye solid (data by Yin et al. 
[TH). The intercept of AP/T^ vs. extrapolated from low- 
temperature data points is in agreement with a glassy con- 
tribution from TLS on the order of 100 ppm. The arrow 
indicates the onset of deviation from the Debye solid. The 
large scatter of data points at lowest temperatures is due to 
the subtraction of a constant pressure Pq. Dashed lines are 
guides to the eyes. 



00 for T < Tg, which describes non-Arrhenius 



the maximum size of Sea is independent of frequency and 
is controlled by eg and qq. 

The results of the acousto-optical model for the dielec- 
tric function are shown in Fig. [TJ We obtain excellent 
agreement with the measurements by Yin et al. 11[ for an 
applied alternating voltage at 500 Hz. Additionally, we 
predict an accompanying dissipation peak in the dielec- 
tric function or a phase lag — arg{e) between the real 
and imaginary part of e{ui,T). Experimental confirma- 
tion of the dissipation peak will constitute an important 
test of the glassy nature and the acousto-optical coupling 
of defects in solid helium. In our calculation, we assumed 
the Vogel-Fulcher-Tammann (VFT) form for the glassy 
relaxation time: t{T) 
r{T) 

thermal activation processes for Tg 7^ 0. Here tq is the 
attempt time and A is the mean activation energy of a 
TLS. In our fit, we are not biasing the ideal glass tran- 
sition temperature Tg to be positive. We find a negative 
Tg, which means no real phase transition occurs at finite 
temperatures, but rather a crossover. This can be due to 
inherent quantum fluctuations originating from the large 
zero-point motion of helium atoms, which suppress the 
onset of a lower phase transition. 

A discussion of the glass model in the context of the 
dielectric function requires a justification for the mag- 
nitude of the elasto-optical coupling to strain. The exis- 
tence of local strain in solid '*He is inferred from measured 
pressure gradients. We obtain a rough estimate, for ex- 
ample, from the mass flow measurements in bulk solid 
''He [13'], where a pressure difference of AP^ ^ 0.1 bar 
across a centimeter-sized pressure cell is reported. The 



estimated local strain, with a bulk modulus B = 320 
bar, is accordingly APl/B = 3 x 10"^. This IS consis- 
tent with the value we used for the fit in Fig. 1, namely 
eo = 8.88 X 10~*. In fact the local strain may have been 
even larger in the experiments by Yin due to the more 
complicated geometry of the measurement apparatus. 

Since the acousto-optical coupling model depends on 
the presence of defects, an important question to ask is 
what is the quality of the helium crystals. For that reason 
we are re-analyzing pressure measurements by Yin [Tlj . 
The glassy contribution of the TLS to the bulk pressure 
through P{T) = Po + agiT'^ + aphT'^, where Pq is the crys- 
tal pressure of bulk solid at T = 0, while coefficients Ugi 
and Qph correspond to the contributions of the TLS glass 
and phonons, respectively. In Fig. [5] we show the temper- 
ature dependence of the deviation AP = P{T) — Pq. A 
large glassy contribution can be seen at low temperatures 
similar to reports by others fl^. The P{T) curve in Yin's 
measurement clearly deviates from a purely Debye lattice 
behavior at around T = 0.4 K with a large positive inter- 
cept corresponding to the order of 100 ppm of TLS 
This number is roughly five times larger than the most 
disordered sample in Lin's IJ] specific heat experiments 
on ultrapure "^He with less than 1 ppb of '^He impurities 
grown over four hours using the blocked capillary method 
[J, [l4| . Clearly crystals grown by Yin and coworkers are 
disordered with sufficiently many TLS defects to support 
centers of local strain fields. 

So far we assumed that both local and global stress 
are constant at low temperatures. From Fig. [5] we can 
read off that the global pressure change between 300 mK 
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and 40 mK is less than APt = 0.18 mbar. This is more 
than three orders of magnitude smaUer than the local 
stress ctl = 8.88 x 10~^ x 320 bar - 250 mbar inferred 
from the dilatory strain cq used in the fit, as well as 
the static pressure difference APl = 0.1 bar measured 
at two pressure gauges in the mass flow experiment |13|. 



Taking all these observations together, we find that the 
change of the dielectric function based on global density 
changes in the Clausius-Mossotti equation is negligible. 
The corresponding density change is Ap/p = APt/B < 
10~^, leading to a change in the dielectric function of 
(5e w (e - I) Ap/p = 0.065APt/S < 10~^ which is more 
than two orders of magnitude too small to account for 
the observed effect ^ 10~^. 

In conclusion, we have shown that the arrested glass 
dynamics causes low-temperature anomalies in strained 
solid ^He through the acousto-optical coupling. We ar- 
gued that the temperature behavior of the dielectric func- 
tion is coupled to local strain fields near crystal defects. 
Thus it records the glassy dynamics and concomitant 
freeze-out of TLS excitations, which also leads to a stiff- 
ening of solid with decreasing temperature. This effect 
is not captured by the standard Clausius-Mossotti re- 
lation, which attributes dielectric function changes to a 
change in mass density or polarizability of the helium 
atom. An important consequence of the phenomenologi- 
cal glass susceptibility is the decrease of dielectric func- 
tion at high temperatures, accompanied by a broad dis- 
sipation peak that can be measured by imaginary part 
of dielectric function. We hypothesize that the coopera- 
tive motion of atoms forming the TLS along dislocation 
segments is the relevant process contributing to the re- 
ported anomaly. In our model, both the change in e{uj, T) 
and the dissipation peak are independent of applied fre- 
quency. Since the glass parameter of the back action 
depends on the concentration of the TLS, we predict that 
the change in dielectric function will be larger in quench 
cooled or shear stressed samples. 
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